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1. Introduction 

One of the least documented aspects of the massive star formation feedback is the role they 
might play in accelerating cosmic rays (CRs). Indeed, the primary sites for CR acceleration mech¬ 
anisms have long been considered to be shocks associated to rather young supernova remnants 
(SNRs, e.g. [8]). However, [24] was the first to seek a correlation between star formation (traced 
by OB associations) and y-ray emission as detected by the Cos B satellite. In this scenario, the OB 
association served the purpose of low-energy particle injection, and the neighbouring, expanding 
supernova remnant (SNR) provided the necessary energy to accelerate them. In turn, [39] inves¬ 
tigated the acceleration of cosmic-rays at the terminal shocks of OB star winds, [41] studied the 
subsequent synchrotron emission due to particle acceleration in hot stars environments, and [29] 
showed evidence for non thermal phenomena in Wolf-Rayet stars at X- and y-ray frequencies. 
Finally, [33] demonstrated a significant evidence for an association of y-ray sources detected by 
EGRET with OB-star forming regions (SERs). High-energy telescopes were then turned towards 
massive SERs, leading to the detection of TeV radiation with HEGRA in the vicinity of the Cygnus 
OB2 association region ([3]), and above 380 GeV with HESS in the young stellar cluster Wester- 
lund 2 ([4]). More recently, the detection of TeV emission with Fermi was reported in the massive 
SER W43 by [21], and in the Cygnus X superbubble by [2], more or less simultaneously to the 
study of the association of Fermi sources with massive young Galactic objects by [27]. Erom the 
modelling point of view, in the massive SER IRAS 16547-4247, the observed spectral index was 
found to be consistent with the existence of an uncooled population of relativistic electrons pro¬ 
duced by diffusive shock acceleration at a strong non-relativistic shock. This triggered the study of 
[6] and [10], who concluded to the possible detection in this region of signatures of accelerated CRs 
by the Cherenkov arrays such as HESS-II at very high energies {Ey > 100 GeV), in an affordable 
amount of time. An overall conclusion is that it is possible that all favourable conditions for CR 
acceleration be co-existing in massive star-forming regions, where strong magnetic field, a high 
degree of ionization, and mechanical energy in fhe form of shocks are simulfaneously presenf, pre¬ 
cisely owing fo fhe effecls of shocks and HII region. In fhis confribufion, we will only summarise 
fhe resulfs of our sfudy of fhe irradiafed shocks in fhe W28 A2 massive sfar forming region, based 
on observafions from fhe APEX, IRAM 30m, and Herschel felescopes, and presented in Gusdorf 
el ah, recenfly submitted lo A&A. 

2. The W28 A2 complex star-forming region 

W28 A2 is a high-mass sfar forming sife of complex morphology. Ils sfruclure and intersfellar 
conlenl is besf described by fhe mosl recenl observafions (e.g. [20]), which benefited from better 
sensilivily and/or spatial resolulion lhan Iheir older counlerparls (e.g. [19]). A skefch of fhe region 
can be found in Eigure 1, indicaling fhe main componenls of fhe region: fhe cenlral sfar or sfar clus¬ 
ter, fhe sun'ounding HII region and associaled dense cocoon of gas and dusl, oulflows revealed in 
various Iracers, maser spofs, X- and y-ray emission peaks. We adopted fhe dislance delerminalion 
of 1.28j(QQg kpc by [26], oblained Ihr'ough aslromelric observafions of fhe 22 GHz water maser 
wilh VERA. Al fhe cenlre of fhe region, an ionizing sfar was direclly imaged by [14]. Given fhe 
larger dislance Ihese aulhors associaled fo fhis source, [26] revised ifs classificalion in favour of a 
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Figure 1: A symbolic sketch of the W28 A2 region. The grey circle marks the beam of our analysis. The 
red star is the central, ionizing, so-called Feldt’s star. It is surrounded by an orange circle that delimitates the 
extent of the ultra-compact HII region, itself surrounded by a dusty and molecular shell shown in the form of 
a green circle. The blue crosses indicate the positions of five continuum sources. The purple star shows the 
centre of a Bry outflow. The yellow triangles indicate the H 2 knots attributable to outflows apices. The two 
corresponding outflows are shown in the form of dark blue ellipses. The spectacular outflow observed by 
[19] is shown as light blue ellipses. Finally, maser spots are indicated; NH 3 (green squares). Class I CH 3 OH 
(pink circles). Not shown are the H 2 O and OH ones. Adapted from [20], shown in Gusdorf et al., subm. 


08-08.5 ZAMS. Around this ionizing star, [20] report the detection of five dust-emission sources, 
that could be protostars or protostellar cores at various stages of evolution. The central star has 
carved a roughly circular ultra-compact HII region (e.g., [25], [42]) surrounded by a dense cocoon 
of dust and molecular gas (e.g., [36]). The central region resembling a cluster of sources, drives a 
gigantic East-West outflow first reported by [19]. Other outflows have been identified in the region: 
two small, nearly North-South ones were detected through interferometric observations of CO and 
SiO lines (e.g. [40, 20, 37]), and their apices could coincide with the observations of H 2 knots at 
2.12/im by [30]. Finally, a small outflow has been detected west of the ionizing, central star, in the 
Bry line by [30]. The presence of numerous maser spots has also been reported in the region in var¬ 
ious species, in OH, H 2 O, CH 3 OH, and NH 3 (see e.g. [20] and references therein). [15] and [35] 
used their OH maser observations to constrain the local strength of the magnetic field to around 
2 mG. This value was also measured on a region of large spatial extent (>10000 A.U.) by the ap¬ 
plication of the Chandrasekhar-Fermi method to the high angular resolution SMA observations of 
the polarised dust continuum emission at 870 /^m by [38]. 

W28 A2 has also been widely observed in the high-energy regime owing to its proximity with 
the W28 SNR, and detected independently from it. Although not seen in the X-ray ROSAT and 
ASCA detector ranges (respectively 0.5-2.4 and 0.6-10 keV) by [32], the HII region was comprised 
in the peak of emission detected by XMM-Newton, and was very clearly detected by Chandra, 
according to the preliminary results presented by [34]. Furthermore, its location seems to match 
that of one very high energy (F > 0.1 TeV) y-ray emission peak described in [5]: HESS J1800- 
240B, suggesting that CR acceleration might be taking place. This view is supported by the most 
recent y-ray observations of the region by AGIFE (E > 400 MeV, [17]), and by the Fermi FAT 
(from 0.2 to 100 GeV, [1]), which both confirmed prominent emission on this location. The Fermi 
and HESS combined y-ray spectra are accounted for by [16], [22] and [18], who argued that the 
strong energetic y-ray emission is caused by the collisions from energetic protons (escaping from 
the SNR shock front) on the W28 A2 molecular cloud, due to the close distance between the high- 
mass SFR and the W28 SNR. However, the latest study recognized the lack of models to evaluate 
the possible in situ acceleration of CRs within the HII region or the outflows associated to this SFR. 


3 





Cosmic rays acceleration in the shocks of the W28 A2 massive star-forming region ? 


Antoine Gusdorf 




Figure 2: Left panel: a sample of our CO lines observations; (6-5) and (7-6) (green and blue lines, APEX 
telescope), (16-15) (red line, Herschel telescope). Central panel: the CO rotational diagram, obtained 
from Herschel observations (red points). Right panel: atomic carbon transitions observed with the APEX 
telescope: ^Pi-^Pq and ^P 2 -^Pi (black and red lines). All data were obtained in the central position of 
W28 A2 (see Figure 1), and extracted from Gusdorf et al., subm. 

3. Evidence for irradiated shocks in W28 A2 

We mapped the emission of the rotational transitions CO with /up = 3, 4, 6, 7 with the APEX 
telescope, performed pointed observations of the /up = 6, 7, 16 with the HIFI receiver onboard 
the Herschel telescope, and obtained a footprint of the full ‘Spectral Line Energy Distribution’ 
(SLED) accessible to the PACS receiver onboard the Herschel telescope. We also observed the 
*^CO isotopologue’s transitions /up = 3 and 6 with the APEX telescope, while the '^CO (14-13) 
line was also detected with PACS. A comparison between the /up = 6, 7 (from APEX) and 16 
(from HlFl) spectra in the central region, at the same spatial resolution, is given in Figure 2. The 
presence of CO shocked gas can be seen in these spectra through the line wings extending up to 
±80 km s^^ towards the blue-and red-shifted directions. From these spectra, we also infeiTed that 
approximately one third of the gas is comprised in each kinematical component: blue-shifted (-100 
to 4 km s^'), ambient (4 to 16.5 km s^^), and red-shifted (16.5 to 100 km s^'). From our maps, 
we also concluded that emission in all CO lines was filling the beam of our observations (12.5", or 
0.08 pc); we hence adopted a filling factor of 1 for fhis emission. However, our observafions from 
fhe ground also showed fhaf all lines af leasf up fo /up = 6 were opfically fhick, even in fhe wings. 
Our PACS observafions, on fhe ofher hand, demonsfrafed fhaf fhe CO (14-13) line is opfically fhin, 
which we exfrapolafed fo fhe emission from fhe higher-lying transitions. Using those lines, we 
could then resort to a rotational diagram approach (shown in Figure 2), and derived the excitation 
conditions for CO: the column density in each kinematical component was found to be of the order 
of lO'^ cm^^, with a minimum kinetic temperature of the order of 270 K. 

Similarly, we mapped the ^Pi-^Pq and ^P 2 -^Pi line emission of [Cl] with APEX, with their 
line wings also broadly filling a beam of 12.5". However, fhese wings were found fo be much less 
brighf fhan CO ones, as can be seen in Figure 2. We performed an LVG analysis fo infer a column 
density of (1.5-3) xlO'^ cm^^ in each blue- and red-shiffed componenfs, wifh respective kinetic 
femperafures of 140-280 and 55-70 K. These laffer values are subjecf to fhe calibrafion uncertainly 
on fhe (2-1) line wing emission. 

Additionally, we mapped fhe [Cll] line (^P 3 / 2 -^Pi/ 2 ) 1900 GHz wifh fhe Herschel fele- 

scope. Again fhe emission was found fo roughly fill a beam of 12.5". The velocify-resolved profile 
exfracfed in fhe cenfral position shows wings exlending up to ±50 km s^\ and conslilules unam¬ 
biguous evidence fhaf fhe shocks in fhe region are irradiated by fhe UV field of fhe prolo-slar(s) 
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Figure 3: Left panel: the spectrum (black line and grey histograms). Right panel: the CH+ spectra; 1-0 
(red line), and 2-1 (black line and grey histograms). All [CII] and CH+ spectra are overlaid with the CO 
(16-15) and [Cl] ^Pi-^Po ones from Figure 2 (blue and pink lines); they were obtained with the Herschel 
telescope in the central position of W28 A2 (see Figure 1), and extracted from Gusdorf et al., subm. 


(the corresponding spectrum can be seen in Figure 3). LTE conditions, optical thinness, and an 
excitation temperature in the range 50-250 K, consistent with the coincidence between [CII] and 
CO (16-15) line wings were assumed. Our analysis then yielded a column density value of (0.5- 
1.2)xl0'* cm^^ in the blue-shifted component, and a lower limit of (1.5-3.5)x 10^’ cm^^ in the 
red-shifted component. This result quantitatively supports the view of a significant role of the 
UV radiation of the forming star, since it makes C+ the dominant form of carbon in the observed 
shocks. Assuming ?i(C)/?i(H)=[?i(CO)-i-?i(CI)-i-?i(C+)]/n(H)=1.4x 10^“^, we found an H equivalent 
column density of (l.l-2.2)xl0^^ cm^^ and (0.6-1.l)x 10^^ cm^^ in the blue- and red-shifted 
kinematical components. This allowed us to estimate a lower limit to the ionization fraction in the 
observed region, given by [n(C+)/n(H)], and of the order of (0.2-1)x 10^"*. This range is only a 
lower limit, as it is based on the assumption that C+ is the dominant charge earner in the medium. 

This also allowed us to infer the energetic parameters associated to the part of the outflows 
caught in the beam of our observations, following the method of, e.g. [9]. The total outflow mass 
caught in the beam of our observations was found to be 0.6-1.2 Mq. The assumed outflow velocity 
in each blue- and red-shifted directions and associated to most of this mass was 50 km s^*, as 
observed for C+, corresponding to a dynamical age of 760 years. This resulted in a total (blue- 
and red-shifted wings) mass-loss rate of (0.8-1.2)x 10^^ Mq yr^\ a total momentum of (32- 
62)Mq km s^^, or a mechanical force of (4-8)x IO^^Mq km s^^ yr^^ Additionnally, we measured 
that a kinetic energy larger than (1.6-3.1) x lO'*^ erg is dissipated in these outflows, corresponding 
to more than 175-340 Lq. These values are rather typical of bipolar outflow activity driven by an 
0-type star or cluster containing one ([23]). 

Finally, the in'adiated nature of the shocks associated to the bipolar outflows was confirmed 
by the observation of the (1-0) and (2-1) lines of CH+ in single-pointing mode with the HIFI re¬ 
ceiver onboard Herschel telescope (shown in Figure 3). Both line profiles exhibit wings extending 
in the red-shifted direction to 80 km s^^ Because of the different beam sizes of these two transi¬ 
tions, we resorted to a simple LTE and optically thin analysis to yield the column density of (4.4- 
12.7)x 10*^ cm^^ in the red-shifted component (corresponding to a fraction of (0.4-2.l)x 10^^ of 
H equivalent). This result is based on an excitation temperature of 50-250 K compatible with the 
assumption that CH+ is present in the same gas component as CO, Cl and C+, as suggested by the 
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similarity of the wing emission of all their observed transitions. This high value for dense gas-CH+ 
is in the range of that derived in a similarly massive star forming environment, DR21 by [13], who 
interpreted it by means of ‘irradiated’ shock models. 

The presence of shocks hence large mechanical energy budget on a rather large scale, where 
the ionization fraction is not negligible and where the magnetic field is intense ([38]) led us to start 
working on possible in situ scenarios of CR acceleration in this region. 

4. In situ cosmic rays acceleration in W28 A2 

As we have seen in Section 2, the W28 A2 complex is coincident with gamma-ray sources 
detected by Fermi and HESS. One possible explanation of the gamma-ray signal is that cosmic 
rays accelerated at the forward shock of the nearby SNR W28 escape and then interact with the 
dense material associated with W28 A2 [16, 22, 18]. Here, we rather study the possibility that 
in-situ acceleration of energetic particles are associated with a part of the observed gamma-ray 
radiation. Additionally, we only worked on the evaluation of the maximum particle energies, with 
views to determine if it could be high enough to at least partly account for the HESS and Fermi 
detections. The work presented in Gusdorf et ah, submitted, added to the literature, and specially 
to those of [26], and [38], allowed us to better understand the complex structure of the W28 A2 
region. In particular, we could highlight three regions where particle acceleration could take place: 

• (i) in the irradiated shocks associated to the large-scale bipolar outflows, characterized by 
a shock velocity f7sh=50 km s^*, age tage=760 years and a radius of 16000 AU with an 
ionization fraction of ~ 10^^ (all inferred from our carbon species observations), an ambient 
density uh = 10^ cm^^ and temperature r=100 K (typical of shocked gas densities and 
inferred from our observations), a magnetic field B=2 mG ([38]); 

• (ii) in the HII region, assuming the outflows are launched in this region close to the central 
star, with the following parameters: shock velocity f/sh=50 km s^^ and age tage=600 years 
(both inferred from our carbon species observations), an ambient density uh = 5 x 10^ cm^^, 
a magnetic field B=1 mG ([38]), a temperature of r=10000 K typical of HII regions, and a 
radius of 3300 AU ([26]); 

• (iii) in the central cluster, where the supersonic winds of the massive stars can interact and 

accelerate particles, then injected in the HII region. This scenario is more speculative but 
might be the most efficient. The central star is of O spectral type, and powers winds with 
typical terminal velocities ~ 2000 — 3000 km s^^. In the following, we use the inner radius 
of the HII region estimated by [26] as the typical size: 1200 AU. We also assume an age 
of 10^ yrs typical of forming stars surrounded by an ultra-compact HII region. Based on 
a mass-loss rate of yr^' (more than two orders of magnitude less than the value 

associated to the bipolar outflow activity and inferred from our observations), the resulting 
kinetic energy thus deposited by the stellar wind is of the order of 10^*^ ergs. 

Eor the scenario (i), we found that CR acceleration was not very likely in the outflow region. It 
is indeed weakly ionized with a typical ionization fraction ~ 10^“*. Neutrals can strongly limit the 
maximum particle energies through ion-neutral damping of the resonantly interacting waves ([12]). 
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We found that ion-neutral collisions are so strong that they can even prevent particle acceleration 
to enter the relativistic regime. Much faster shocks are needed to generate any y-ray signal from 
these regions. Regarding scenario (hi), inspired by [28, 11], we found that the maximum particle 
energy is limited by escape losses. We also measured that if a small fraction of the kinetic energy 
driven by the stellar winds was imparted into energetic particles, it could be sufficient to account 
for the energy necessary to explain gamma-ray emissions. It is also possible that such a star is 
surrounded by several less massive ones. In that case, energetic particles can be re-accelerated by 
multiple shock interactions. We can speculate that an acceleration by the central star/cluster is a 
possible mechanism of injection of high energy particles in the ultra compact HII region, but this 
has to be constraint further by additional observations and modelling efforts. 

Finally, in the HII region (scenario (ii)), the sonic and Alfvenic shock Mach numbers can 
be estimated: Mg ~ 5.4 and ~ 8.1 respectively, and show that the shocks are strong and super- 
Alfvenic. The region being fully ionized the maximum particle energy is limited either by the age of 
the shock or by losses. Particles can be lost for the acceleration process either because they escape 
upstream once their diffusive length is of the order of a fraction of the shock radius or as they loose 
their energy radiatively [7]. The radiative losses depend on the type of particle: electrons loose their 
energy either by synchrotron or Inverse Compton losses whereas protons or heavier nuclei loose 
their energy through inelastic collisions with the surrounding matter. Following the work of [31], 
we derived the maximum particle energies depending on the dominant energy limiting process. 
Regarding hadrons, we found that the main limiting factor of the hadron energy in parallel shocks 
are geometrical losses upstream the shock with typical maximum energies Fmax.h ~ Z x (5 — 135) 
GeV. As for the case of electrons, the radiative losses strongly limit the maximum electron energies 
to values Fmax.e ~ 5 GeV. It is not likely to obtain a TeV signal with these particle energies as 
Inverse-Compton scattering of more energetic photons rapidly fall in the Klein-Nishina regime. 
We hence overall concluded that unless invoking local magnetic field amplification, faster shocks 
or/and quasi-perpendicular shock configuration, the maximum hadron energies do not reach the 
multi TeV energy domain and can hardly explain the HESS signal. However the maximum particle 
energies are compatible with the Fermi waveband. The low maximum electron energies and the 
fact that the highest soft photon energies are in the UV band does not favour a leptonic scenario. 
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